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Abstract 

Average multiplicities and transverse momenta of hadrons are cal- 
culated using a microcanonical hadronization description for a cluster 
of given total energy and volume. As a function of the total energy, 
we determine the critical volume above which the microcanonical de- 
scription coincides with the canonical one, and compare the results 
with those obtained using one of standard canonical models. We show 
that the critical volume depends on the energy and the mass of the 
hadrons. For heavy particles, volumes above 50 /m^ are needed, even 
more thanlOO/m'^ if one considers transverse momenta. Thus the pre- 
diction of heavy hadron multiplicities in pp, Kp, and e+e^ reactions re- 
quires a microcanonical approach, whereas for heavy ion reactions a 
canonical calculation is valid. We conclude by showing the importance 
of the feeding for the observed hadron multiplicities. 



1 Introduction 

The statistical description of proton-proton and heavy-ion reactions has al- 
ready a long tradition. It was Hagedorn (Tl El, who showed in a sequence 
of papers that many aspects of these reactions are close to that one expects 
assuming that the transition matrix element is constant and therefore the 
distribution of final state particles is completely determined by phase space. 

More recently the statistical interpretation of nuclear reactions regained 
interest, after it had been demonstrated that in heavy ion reactions at 
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CERN-SPS amd RHIC energies the multiplicities of a multitude of non strange 
hadrons is in remarkable agreement with the assumption that all particles 
are created in thermal equilibrium at a temperature which is close to that 
expected from lattice calculations for the transition of a hadron gaz towards 
a quark-gluon plasma [3-5]. Also strange hadrons fit into this picture if one 
includes a penalty factor for each strange quark which is contained in a 
hadron. 

Later this approach has been successfully extended towards AGS and 
SIS energies 0, as well as towards pp and e+e reactions (71 |8l. Whereas 
the former are large systems with a large reaction volume, the latter ones 
yield only a small hadron multiplicity and require volumes of the order of 25 
fm^. For such small systems, it is not evident that a canonical or a grand 
canonical description is justified. 

Therefore it is worthwhile to check whether for these small volumes and 
particle multiplicities a microcanonical description still coincides with the 
canonical one. It is the purpose of this article to investigate this question 
employing for the first time a numerical realization of the algorithm to cal- 
culate the microcanonical phase space which has been presented by us in 
ref. EI. For the present investigation, we limit the number of hadrons to 54, 
which include pseudoscalar and vector mesons (octet and singlet) as well as 
the octet and decouplet of baryons and antibaryons. Strange particles are 
produced according to phase space, i.e. without applying any suppression 
factor. The results are compared to that of two canonical calculations using 
the approach of Becattini et al. El |7l but without strangeness suppres- 
sion: in the first calculation the number of hadrons which can be produced 
is limited to the same 54 species allowed in our microcanonical treatment, 
in the second calculation the standard set of hadrons |4, 7, 8| is included. 
In principle we can include more hadrons in the microcanonical ensemble. 
This makes a detailed comparison more difficult, because the less known 
decay channels of the additional hadrons have to agree. 

2 Microcanonical Calculation 

Following the general philosophy of statistical approaches to hadron pro- 
duction, we suppose that the result of a high energy collision can be con- 
sidered as a distribution of "clusters", "droplets", or "fireballs", which move 
relative to each other. Here, we are only interested in An particle yields and 
average transverse momenta, and therefore collective longitudinal motion 
needs not to be considered, and the distribution of clusters may be identi- 
fied with one single "equivalent cluster", being characterized by its volume 
V (the sum of individual proper volumes), its energy E (the sum of all the 
cluster masses) , and the net fiavor content Q = {Nu — Nu, Nd — N^, Ng — Ns). 
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The basic assumption is that a cluster, characterized by V, E, and Q, de- 
cays "statistically" according to phase space. More precisely, the probability 
of a cluster to hadronize into a configuration K = {hi, . . . ,hn} of hadrons hi is 
given by the microcanonical partition function Q{K) of an ideal, relativistic 
gas of the n hadrons hi ||9il, 

= 72^ n / Ud'p^HE - ^Si) 6m) 5q,e., 

^ ^ i=l aeS 1=1 

with Ei = A/mf+~pf being the energy, and pi the 3-momentum of particle 
i. The term ^Q Sg^ ensures flavour conservation; is the flavour vector of 
hadron i. The symbol S represents the set of hadron species considered: we 
take S to contain the pseudoscalar and vector mesons (tt, K, r], r]', p, K*, u, (p) 
and the lowest spin-i and spin-| baryons (A^, A, S, S, A, S*, S*, i7) and the 
corresponding antibaryons. Ua is the number of hadrons of species a, and 
Qi is the degeneracy of particle i. 

We are going to employ Monte Carlo techniques, so we have to generate 
randomly configurations K according to the probability distribution Vt{K). 
We need a method in particular for intermediate size droplets, covering 
droplet masses from few GeV up to 100 or 1000 GeV. So the method should 
work for particle numbers n = \K\ between 2 and 10^, which means, we have 
to deal with a huge configuration space. Such problems are well known 
in statistical physics, and the method at hand is to construct a Markov 
process. So for a given cluster with mass E, volume V, and fiavor Q, we 
start from some arbitrary initial configuration Kq, and generate a sequence 
Ko,Ki,..., Kj^^, with Jeq being sufficiently large to have reached equilibrium 
(which is defined to be the steady state of the Markov process) . If we repeat 
this procedure many times, getting configurations KjI^^,Kj'^^, . . . , these con- 
figurations are distributed as Q{K). We need a transition probability p such 
that it leads to an equilibrium distribution Q{K), with the initial transient 
Jeq being as small as possible. Such an algorithm has been realized for the 
first time in |91 . 

The problem is solved in several steps. One first writes the phase space 
integral as 

/n 
l[£p,6iE-i:E,)6m) (1) 
i=i 

/n n n 

n dPiWpi - E ^(pi' • • • ' P")' 
i=l i=l i=l 

with Pi = \pi\, and with the "random walk function" W given as 

^ „ n n 

W{p,,...,p^)■.= J—^ \[d^,5{Y,P^e^). (2) 

^ ' ^ 1=1 i=l 
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with ii = Pi/lpil- The name "random walk function" is due to the fact that 
W represents the probability to return back to the origin after n "random 
walks" piii with given step sizes pi. In |9|, many details can be found about 
an efficient calculation of W for any n (big or small) . 

The next step amounts to getting rid of the energy delta function. A 
variable transformation gives 



(f){E,mi, . . . ,mn) (3) 
dri... / dr„_i mi,..., m„;ri,...,r„_i), 



JO 



with 

?/'(E,mi, . . . ,m„;ri, . . . ,r„_i) (4) 



UP^e,W{p„...,Pr.). 



1=1 



The symbol T denotes the total kinetic energy E — Hrrii, and the absolute 
values of the momenta are expressed in terms of the as 



Pi = \/U{ti + 2mi) 

ti = T(Xi — Xj„i), Xq = 

Xi = Xj+iA/ri, Xn = 1. (5) 

In principle one may use Monte Carlo techniques to calculate the inte- 
gral, but this is very time consuming. A more elegant method amounts to 
generalizing the hadronic final state by considering not only hadron species, 
but also their momenta. So we use the generalized configurations 

G = {/ii,...,/i„;ri,...,r„_i}, (6) 

where the are related to the momenta pi via eq. O. The weight of such a 
configuration is 

^^^^ " n n "'I' • • • ' • • • ' ^"-i)- 

^ ' i=l a£S °'' 

This expression is well suited to generate configurations G according to 
Q{G), by constructing Markov chains, for details see IQJ. 

Flavor conservation is trivial to take into account, by considering only 
propositions in the Markov chain construction which conserve the total fla- 
vor Q. 

Our algorithm provides a fast method to generate hadron configurations, 
characterized by the number of hadrons, their type, and their momenta. 
In this sense we have a real "event generator" of statistically generated 
hadrons. 
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3 Results 



3. 1 Multiplicities 

In order to see best how the results of a microcanonical calculation ap- 
proach those of a canonical one, we investigate the average particle density 
Pa = na/V, where is the average particle multiplicity of hadron species 
a, and V is the volume. In the microcanoncial approach, the average parti- 
cle density is in general a function of the two variables E and V . For large 
volumes, however, p should only depend on the ratio e = E/V, and this rep- 
resents the limit where microcanonical and canonical description should 
coincide. 

We therefore calculate the density pa as a function of the energy density e, 
for different volumes. We expect these curves to converge for large volumes. 

Figs. 1-3 show the average particle density pa as a function of the energy 
density e for three different volumes V= 12.5 (line), 25 (dashed line) and 50 
(dotted) fm^ The baryon density is fixed at 0.08 haryons/ fm? and the total 
charge is 2 for reasons which will be explained later. The particle densities 
for V = 12.5 /m^ and V = 25 fm^ differ at low energy densities for heavy 
mesons as well as for many of the baryons. The deviation increases - as 
expected - either with increasing particle mass or with increasing strange 
quark content, because the strange quarks have to be compensated by other 
particles in order to obtain the quantum numbers of the two incoming pro- 
tons. 

We see as well that all mesons and almost all of the baryons have already 
arrived at the canonical limit for a volume as small as V= 25 fm?. The only 
exceptions are the heavy triple strange barons Vt and Vt at low energy den- 
sities, where the total energy of the disintegrating system is (at V= 25/m^) 
only a couple of times larger than the mass of the Vt ,0 pair. 

From this observation we can conclude that the canonical limit is ob- 
tained if the total energy of the system is about 10 times the energy of the 
sum of the particles under consideration and of that particles necessary to 
compensate the deviation of the quantum numbers of the considered parti- 
cle from that of the total cluster. 

With the exception of the proton the density of all particles increases 
with increasing energy density. For the proton this is not the case because 
at very low energy densities the baryon cannot get rid of its charge, whereas 
at higher energies a produced meson can carry this charge leaving behind 
a neutral baryon. 

In Figs. 1-3 we have marked as a dot the results of a canonical calcula- 
tion with the parameters fitted to describe the particle yields observed in pp 
collisions at 27.4 GeV |4l 0. The parameters used in this calculations are 
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Figure 1 : Density of particles as a function of the energy density of a cluster 
for three different volumes: 12.5 fm^ (full line), 25 fm^ (dashed line) and 50 
fm^ (dotted line) for an initial baryon density of 0.08 baryons / fm^ and a 
total charge of 2 using a microcanonical phase space calculation. The dots 
present the result of a canonical calculation provided by Becattini. In both 
cases the number of hadrons is limited to 54 and strange particles are not 
suppressed. 
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Figure 2: Same as previous figure, but for additional hadrons 



V = 25.5 fm^ and T = 162 MeK. In contradistinction to the calculation with 
these parameters which is presented in ref. |7|, here the strange particles 
are not suppressed by a 7^ factor and only the 54 hadron species mentioned 
above are produced. Therefore this canonical calculation can directly be 
compared with our microcanonical approach. The average energy of the 
clusters obtained in this calculation is 8.74 GeV, resulting in an average 
energy density of e = 0.342 GeV/fm^. We observe that for all light hadrons 
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Figure 3: Same as previous figure, but for additional hadrons 



the agreement is very nice, verifying that indeed at a volume of y = 25 /m^ 
the multiplicity has arrived at its canonical limit. For the heavy particles, 
such as and Vt, volumes above 50 frr? are needed. 

Another interesting question is how the observed particle yield is influ- 
enced by feeding. It is addressed in flgs. 4-6 which show for the 54 hadrons 
two densities: that of those particles which are produced directly and that 
present after electromagnetic and strong decays. Both microcanonical cal- 
culations are compared with the corresponding canonical results. We see 
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Figure 4: Density of particles as a function of the energy density of a cluster 
with a volume of 50 /m^ in the microcanonical calculation (lines) in compar- 
ison with a canonical calculation at a energy density of .342 GeV/ fiv? with 
54 hadrons (solid points) and 250 hadrons (open points), before (solid lines 
and squares) and after (full lines and circles) strong and electromagneticc 
decay . 

that feeding is unimportant for the vector mesons, whereas for the pseu- 
doscalar mesons feeding increases the yield by more than a factor of 2 for 
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Figure 5: Same as previous figure, but for additional hadrons 



the tt's and by more than 50% for the K's. With the exceptions of protons, 
neutrons and A's feeding is much less important in the baryonic sector. The 
strange (anti)baryons in the electromagnetic decay chain of Q and Q show 
some feeding whereas for A's, fi's, x*'s and S*'s feeding is absent. We display 
in figs. 3-6 as well the results of the canonical calculation with the same 
parameters in which the standard set of hadrons is included. For the heavy 
hadrons both canonical results differ little, whereas the standard set yields 
more pions and kaons, the decay products of the hadrons not included in 
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Figure 6: Same as previous figure, but for additional hadrons 



54 particle set. For the pions both canonical results differ by less than 
50\% for all the other hadrons the differences are much less important. 

3.2 Transverse momentum 

In fig. 7, we plot the average transverse momentum < > as a function of 
the hadron mass m. We compare our microcanonical results (points) with 
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Figure 7: < prim) > for different droplets with an energy density of 
0.34:2 GeV / fm^ as compared to the value for the equivalent canonical sys- 
tem (T=162 MeV). On the left hand side we display the results for droplets 
with baryon numbers B=0 and B=2, respectively, the latter with and with- 
out strong decay. One the right hand side we present those for a droplet 
with different volumes and with a baryon number B=0. 

an energy density of e = 0.342 GeV/fm^ to the canonical result, using 

(vrmT/2)(V^)ir5/2(m/r) 
< ^) K^MT) 

El [T0|, with a temperature of T = 162 MeV (solid lines). The < px >values 
obtained in the canonical calculation of Becattini agree with that analytical 
formula. 

On the left hand side, we display < > of hadrons before decay as 
a function of their mass, for a volume of V = 25 fm^, and a total baryon 
number B = (stars) and B = 2 (squares). Whereas the < pr > of n's and 
K's are close to the value one expects for an equivalent canonical system, 
for the heavier particles the deviation from the canonical value increases 
with increasing mass. The reason is easy to understand: If such a heavy 
particle is produced, the available energy is not sufficient to fill phase space 
up to the high momenta. We see as well that the < px > of the Q is 1.5 
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time the value of a tt . Please note that relativity moderates the increase. 
In a nonrelativistic calculations we expect that <pt > increases as y/m. We 
display on the left hand side as well the average transverse momenta after 
the strong electromagnetic interaction decay for the case of B — 2 (cycles). 
The decay modifies considerably the < pt >of the light mesons and baryons. 

On the right hand side, we study the critical volume for a droplet with 
total baryon number (in order to keep the baiyon density constant when 
changing the volume). We see that with increasing volume the < > of 
the heavy baryons approaches only very slowly the canonical value. How- 
ever, even at volumes of 100 fm^ the asymptotic value is not reached. Thus 
kinematical observables are much more suited to observe the deviations 
between canonical and microcanonical systems. 

4 Conclusions 

We have presented the first numerical calculations of the multiplicity distri- 
bution of hadrons assuming that they are distributed according to micro- 
canonical phase space. In this study we have limited ourselves to a set of 
54 hadrons. The results have been compared with the results of a canonical 
system which has the same size and the same average energy, and agree- 
ment has been found. We observe that for the multiplicity the canonical 
limit is obtained, if the energy of the system is around ten times the mass 
of the observed particle plus that of those particles which are necessary to 
compensate its deviation from the quantum numbers of the cluster. These 
results raise doubts whether heavy hadrons in small systems like pp, Kp 
cind e+e~ can be described in a canonical approach but shows as well that 
for a heavy collision a canonical treatment is valid. Both, canonical and 
microcanonical calculations, show that feeding is very important for the 
pseudoscalar mesons as well as for p,n and A but that is of minor or no 
importance for the other particles. We observe, as predicted analytically, an 
increase of the average transverse momentum < pr > with the mass of the 
hadron. However, we also find that a very big volume is required to obtain 
the value of the equivalent canonical system. 
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